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Abstract -We have studied variability of the spectral lines of the OB star S Ori A-the brightest component of 
the S Ori triple system. Forty spectra with signal-to-noise ratios «500-800 and a time resolution of four minutes 
were obtained. We detected variability in the Hell A 4686, Hel A 4713 and H/3 absorption and the CIII A 5696 
emission line profiles. The amplitude of the variability is «(0.5-l)% of the continuum intensity. The dynamical 
wavelet spectrum of the profile variations reveals large-scale components in the interval 25-50 km/s that move 
within the -Vs'mi to V sinz band for the primary star of the system, Aa 1 , with a band crossing time of 4^-5^. 
However, some of the variable features go outside the band, presumably due to either imhomogeneities in the 
stellar wind from S Ori Aa 1 or non-radial pulsations of the weaker components of the system, Aa 2 or Ab. The 
detected variability may by cyclic with a period of rj A h . We suggest that it is associated with non-radial 
pulsations of the primary in the sector mode (I, m) — (2, —2). 



1 Introduction 

Spectral observations of hot stars at UV PQ [2] , op- 
tical [DEI El El HE], and X-ray [3 E] wavelengths 
indicate the presence of structures in the atmo- 
spheres of these stars with various sizes and den- 
sities, and with lifetimes from less than an hour to 
several days. Variations of the line profiles in the 
spectra of OB stars are mainly regular or quasi- 
regular. The spectra of B stars and six O stars 
display regular short-period (3 - 12 h) profile vari- 
ations for numerous lines (in particular, Hel) as- 
sociated with the non-radial pulsations of these 
stars [BJ. 

Weak magnetic fields (several hundred Gauss 
at the stellar surface) may be one factor favor- 
ing the formation of large-scale structures in the 
atmospheres of O stars [S]. Currently, magnetic 
fields have been detected with sufficient reliabil- 
ity only in one O star, 9 1 Ori C [10], and several 
early B subtype stars[llj. Unlike Bstars, the am- 
plitudes of the line profile variations in the spec- 
tra of O stars are relatively small (l%-3%; see, for 
example, [12]), so it it is more appropriate to call 
this line profile microvariability. To detect such 
variability and to clarify its nature, observations 
with high time and spectral resolution and high 
signal-to-noise ratios (> 300) are needed. 

Here, we continue our study of the microvari- 
ability of hot stars begun in [12], when we ana- 
lyzed line- profile variability in the spectra of the 



supergiants 19 Cep, a Cam A Ori A (09.511). 

We now present the results of spectral observa- 
tions of the spectral triple system 5 Ori A, whose 
primary displays the same spectral type as A Ori A 
- 09.5. Section 2 contains general information 
about the 5 Ori A system, and Section 3 describes 
the observations and reduction of the spectra. Sec- 
tion 4 presents the methods used in and the re- 
sults of our analysis of line-profile variability in 
the spectrum of 5 Ori A.. Our interpretation of 
the observations is presented in Section 5, followed 
by a brief conclusion section. 

2 THE 5 Ori A SYSTEM 

5 Ori is a wide visual triple system contain- 
ing three components: A (HD 36486), B (Bd- 
00°.983B), and C (HD 36485), with apparent 
visual magnitudes of 2.23 m , 14.0 m , and 6.85 m , 
respectively. Components Band C are at dis- 
tances of 33" . and 53" , respectively, from the pri- 
mary component, A. Our observations refer to the 
brightest component 5 Ori A. 

5 Ori A (HD 36486, HR 1852) is a phys- 
ical triple system with the primary 5 Ori Aa, 
which is an eclipsing binary with orbital period 
P =5. 73 d [13], and the secondary 5 Ori Ab, with 
an orbital period of 224.5 yrs. Recent detailed 
Doppler tomography studies [14] have made it pos- 
sible to refine the parameters of the .Ori A triple 
system, presented in the Table. 
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Table 1: Parameters of the <5 0ri A system 



Components 



Parameter 


Aa 1 


Aa 2 


Ab 


References 


Distance to the star, pc 


360 


360 


360 


El 

i — i 


Distance to component Aa 1 , Rq 




33 


w 25 000 




Spectral type 


09.511 


B0.5III 


B (early subtype) 


on 


Orbital period 




5.7325 days 


w 200 yrs 


d 


Radius, (Rp>) 


11 


4 




HIE] 

! * 1 — i 


Mass, (M ) 


11.2 


5.6 


w 27 


M 


Luminosity. log(L/ Lp,) 


5.26 


4.08 






Teff, K 


33 000 


27000 




i — » i — —i 


lg Q 


3.4 


3.8 




[TBI fT71 


Contribution to optical radiation 


70 


7 


23 


m 

i i 


in Hel A 6678 line region, (%) 










Tasini (km/s) 


157 ±6 


138 ±16 


w 300 


M 




133 






M 


Voo, km/s 


2000 


1500 




nsii2Di 


Mass loss M, M /yr 


1.1 x 10~ 6 


1.2 x 10~ 7 




Hi [20] 



The brightest component, 5 Ori Aa , is a high- 
mass 09.511 star with an intense stellar wind; the 
star's mass-loss rate is log(M) ~ —6.0 and the 
limiting velocity of the wind is VJ>o ~ 2000 km/s. 
A value a factor of two lower is given in |21j : 
lg(M) w -6.3. According to [22| fl"6]. sa 
2300 /. The rotational velocities of the star de- 
termined by different authors differ appreciably. 
In [14J, a relatively high rotational velocity for 
the primary was obtained, Vsini = 157±6km/s„ 
while Abt et al. [18] found V sini = 133 km/s (see 
the Table). 

The study [16] presents a photospheric radius 
for the primary, Aa 1 R* = 22R.„ twice the value 

= Hi?., given in the Table. This discrepancy 
is probably due to the fact that a larger distance 
to the star was used in [16] than the distance in 
the Table, 450 pc. 

The component 5 Ori Aa 2 is a B0. 5 III star. 
Preliminary studies indicate that the third compo- 
nent of the system, 5 Ori Ab, is an early B-subtype 
star with broad spectral lines, indicating either a 
high rotational velocity (Vsini ~ 300 km/s) or 
binarity of the star p3] • 

Harvin et al. [H] determined the orbital in- 
clination of the 5 Ori Aa close binary to be i ~ 67° 
and estimated the mass of the components to be 
M(Aa 1 ) = 11.2M and M(Aa 2 ) = 5.6M . The 
secondary, 5 Ori Aa 2 , is appreciably weaker than 
the primary (see Table). Both the primary and 
secondary, Aa 1 and Aa 2 , have masses roughly half 
the mass of main-sequence stars with correspond- 
ing positions on the HR diagram [14J.. It was sug- 
gested in [H]. that, in the course of its evolution, 



the system underwent a common- envelope stage 
in which both stars filled their Roche lobe and 
with intense mass loss. After the system had lost 
half its total mass, the distance between the com- 
ponents increased and the mass loss ceased. 

The mass of the third component of the sys- 
tem, Ab, cannot be derived from the radial ve- 
locities. The Table presents the value M{Ab) ~ 
27 Mq from [14], which corresponds to a main- 
sequence star with Mt, i = — 4 m .2. A main- 
sequence star with a mass of 27 Mq should have 
spectral type 08.5V[23j, which contradicts the 
spectral classification in the Table. To explain this 
inconsistency, Harvin et al. [Hj suggest that Ab 
may be a close binary consisting of two main- se- 
quence B0.5V stars with a mass of w 19 Mq. 

3 OBSERVATIONS AND RE- 
DUCTION OF THE SPEC- 
TRA 

Our spectral observations of the 5 Ori A system 
were carried out on January 10/11, 2004 with the 
6m telescope of the Special Astrophysical Obser- 
vatory as part of our program to search for rapid 
line-profile variability in the spectra of early-type 
stars [12]. We used the NES quartz echelle spec- 
trograph [23], which is permanently mounted at 
the Nasmyth focus and equipped with the Upp- 
sala 2048 x 2048 CCD. To increase the spectro- 
graph's limiting magnitude, we used a three-cut 
image cutter [25] . The resulting spectral resolu- 
tion was R pa 60000 and the dispersion 0.033 
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A/pixel. The image size during the observations 
was about 3". The reference spectrum was pro- 
vided by a ThAr lamp. 

The angular separation between the most dis- 
tant components, 5 Ori Aa 1 and Ab, is about 0.3". 
(see Table), which means that all components con- 
tribute to the resulting spectrum of the 5 Ori A 
system. Over a total observation time of ~ 2^50™, 
we obtained 40 spectra of the star with an ex- 
posure time of 180 s. Taking into account the 
CCD readout, the time resolution was 260 s. The 
signal-to-noise ratio per pixel . was 500 in the blue 
(4500 A) and 800 in the red . (6000 A). 

The preliminary reduction of the CCD echelle 
spectra was done in the MIDAS package [T2]. We 
adapted the standard ECHELLE procedures of 
the MIDAS package for data obtained with an im- 
age cutter. The reduction stages included 
1) median filtering and averaging of the bias 
frames with their subsequent subtraction from the 
remaining frames; 2) cleaning of cosmic rays from 
the frames; 3) preparation of the flat-field frame; 
4) determination of the position of the spectral or- 
ders (using the method of Ballester [26J); 5) sub- 
traction of scattered light (to determine the con- 
tribution of scattered light, we distinguished the 
inter-order space in the frames, and then carried 
out a two-dimensional interpolation; this function 
was recorded in individual frames to be subtracted 
from the initial images); 6) extraction of the spec- 
tral orders from the reduced images of the stellar 
spectrum, the flat fields, and the wavelength ref- 
erence spectra; flat-field reduction; 7) wavelength 
calibration using a two- dimensional polynomial 
approximation for the identifications of the lines 
of the reference spectrum in various spectral or- 
ders. 

When the image cutter is used, each order 
of the echelle image is represented by three sub- 
orders (cuts). In a given order, the instrumental 
shift of the upper and lower cuts relative to the 
middle cut was determined via a cross-correlation 
with the reference spectra. The three cuts ob- 
tained in a given order were summed using median 
averaging, taking into account the determined in- 
strumental shifts. 

To study the line profile variability, the pro- 
cessed spectra were normalized to the continuum 
level reconstructed in each echelle order. We 
drew the continuum in spectral orders containing 
narrow spectral lines using the technique devel- 
oped by Shergin et al. [27], in which spectra are 



smoothed with a variable Gaussian filter with a 
25-30 A window, with the same parameters for the 
entire set of spectra. 

In orders containing broad (usually hydrogen) 
spectral lines, the continuum was drawn using the 
following procedure. In all 40 spectra, broad lines 
were cut out in fixed wavelength intervals. To de- 
termine the position of the continuum, we used a 
polynomial approximation for all wavelengths of 
the order, except for the regions of the cut-out 
broad spectral lines. The approximation param- 
eters remained unchanged for all spectra of the 
set. This procedure provides stability and repro- 
ducibility of our determination of the continuum 
in all the spectra with an accuracy of up to sev- 
eral tenths of a percent. This makes it possible to 
reach high accuracy in deriving the differential line 
profiles and detecting profile variability in broad 
lines at levels down to 0.1%. 

4 LINE PROFILE VARIABIL- 
ITY 

4.1 Contribution from Different Com- 
ponents of the System to the Line 
Profiles 

We selected sufficiently strong unblended lines for 
detailed studies of the line-profile variations. The 
selection criterion for absorption lines was that the 
residual intensity be r max > 0.1. These criteria led 
to the selection of the Hell A 4686, Hel A 4713 and 
H^ lines. In addition, we studied the profiles of 
the CIII A 5696 emission line. 

The 5 Ori A system contains three O and B 
stars, each of which could have variable line pro- 
files; the observed profiles result from the contri- 
butions of all stars. It follows from the Table that 
above two-thirds of the flux from the system at the 
wavelength of the Hel A 6678 line comes from the 
primary component, Aa 1 . The contributions from 
Aa 2 and Ab at these wavelengths are substantially 
smaller. 

At other wavelengths, the contributions from 
Aa 2 and Ab may differ from those in the Table. 
To determine to what extent each of the stars in 
the system is responsible for line-profile variability, 
the contribution of each star to the total profiles of 
the analyzed lines must be clarified. In addition, 
taking into account the substantial scatter of the 
Vsini values for Aa 1 (see the Table), reliable v 
sin i estimates for the components of the 5 Ori A 
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system are highly desirable. 




-405 -300 -200 -100 ICS 200 300 400 



V 

Fi gure 1; Hell. 4686 synthetic (bold curve) and observed 
(thin curve) line profiles, plotted in terms of the Doppler 
shift V from the line center. 

To solve these problems, we modeled the com- 
bined profiles for the photospheric Hell A 4686 and 
Hel A 4713 lines in the spectrum of 5 Ori A. We 
assumed the total line profile to be the sum of 
the contributions from the components, since Aa 1 
and Aa 2 are not eclipsed at the above orbital 
phases. We calculated the rotation-broadened 
line profiles using the standard relations (see, for 
example, |28j). The radial velocities of the compo- 
nents for the time of the observations were calcu- 
lated using the ephemeridas of Harvin et al. [T2] . 

During our observations, the orbital phase var- 
ied in the narrow interval 0.848 — 0.851. For such 
small phase variations, the velocities of the com- 
ponents can be considered to be constant over the 
total observation time. The heliocentric velocities 
^(Aa 1 ) = -21km/s, and V(As?) = -102km/s 
were obtained using the ephemeris |14j . 

For convenience in comparison, both the calcu- 
lated and observed line profiles were transformed 
into the system of the center of mass of the close 
binary S Ori Aa. The velocity V T =m 23.0 / [H] 
was used as the zero-point of the velocity scale. 
In this system, the velocity of the third compo- 
nent is V(Ab) =« 27km/s [HE] . 

Figure 1 presents the total profile of the pho- 
tospheric Hell A 4686 line calculated using these 
parameters. The observed profiles are reproduced 



very reasonably. The small discrepancies between 
the observed and calculated profiles near V = 
150 — 250km/s may be due to the contribution 
of the stellar wind from Aa 1 . 

4.2 Variations of the Average Profiles. 
Differential Profiles 

The nightly-average line profiles in O-star spectra 
are often substantially variable [6]. Line-profile 
variations on shorter time intervals are insignifi- 
cant. To make the line-profile variations in the 
spectra of 5 Ori A in time intervals of 30-40 min 
more clearly visible, we divided the profiles into 
four groups with ten profiles in each and calcu- 
lated the average profiles for these groups. 

The average profiles for each group presented 
in Fig. [2] display appreciable variations within an 
hour. These variations are particularly clearly vis- 
ible in the Hel A 4713. In the transition from the 
first group of spectra to the following three groups, 
the line depth decreases by 1%-1.5%. These vari- 
ations are correlated with variations in the groups 
of average profiles for the CIII A 5696 emission line. 
In the transition from the first to the fourth group 
of profiles, the flux in the central region of the 
line increases by ~ l%.The amplitude of the flux 
variations in the H^ and Hell A 4686 lines reaches 
l%-2%. 

In the blue wing of the CIII A 5696 profile, nar- 
row absorption features with depths less than 0.01 
of the continuum intensity can be seen. Similar, 
but less deep, features are also present in the red 
wing of the line. Our analysis shows that the po- 
sitions of all of these features coincide with those 
for atmospheric absorption lines, primarily water 
vapor [29] . 

To reveal variable features in the line profiles, 
we constructed differential profiles by subtracting 
the average profile obtained using all forty avail- 
able spectra of 5 Ori A from the individual profiles. 

Figure presents dynamical differential pro- 
files of the studied lines in the spectrum of 5 Ori. 
Grey shades indicate deviations of individual pro- 
files from those averaged over the entire set of 
spectra. For convenience, the wavelength scale is 
translated into Doppler shifts from the line cen- 
ter. The radial velocity of the center of mass of 
the 5 Ori A system was used as the zero point of 
the wavelength scale (see Section 14. ip . 

We show the differential profiles in the "neg- 
ative", i.e., darker regions in Figj3}x)rrespond to 
regions of the profile that lie above the average 
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A, A A, A 

Figure 2: Average profiles for spectra 1-10 (thin solid curve), 11-20 (dashed curve), 21-30 (bold solid curve), and 31-40 
(dotted curve) for the Hell A 4686 (top left), HeIA4713 (top right), H/3 (bottom left), and CIIIA5696 (bottom right) 
lines. The flux in the line frequencies, F, is given as a fraction of the flux in the continuum adjacent to the lines, F c . 
The arrows indicate the laboratory wavelengths of the lines. 



value (peaks), while lighter regions correspond to 
parts of the profile that lie below the average value 
(valleys). This makes it easier to identify regular 
variations of the profiles in Fig|3J The Hell A 4686, 
HeIA4713, and EL3 line profiles show that the 
peaks move from the violet to the red wing of the 
profile. 

We can see a broad (50-100 km/s), variable 
feature in the vicinity of the Hell A 4686 line. The 



feature appears at a velocity near -100 km/s in 
the first (lower) differential profiles, and is then 
gradually shifted towards the line center. In the 
last (upper) profiles, it is visible in the domain 
of positive velocities. Such behavior of features 
in the differential profiles is typical for non-radial 
pulsations [30]. Virtually no obvious spectrum- 
to-spectrum profile variations near weaker spec- 
tral lines, due to both the short duration of the 
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observations and the low amplitude of the profile 
variations. We used the technique described in the 
following subsection to reveal profile variations in 
such weak lines. 



Analysis of the Spectrum of Time 
Variations of the Differential Spec- 
tra 



4.3 



The method of Fullerton et al. [H] is often used 
to elucidate the presence of weak line profile vari- 
ability. We describe here a substantially modified 
version of this technique, which we used in our 
analysis. Suppose that N spectra of a studied star 
have been obtained. We denote -Fi(A), % = 1, ...N 
to be the flux in the ith spectrum at wavelength 
., normalized to the continuum. Let -Fj(A) be the 
flux at wavelength A averaged over all the observa- 
tions. The dispersion of the random i^(A) is then 
equal to 



N-l 



D(X) = TVS(X) 

n r 

E 9i F i (A)-F i (A) 



(1) 



E 9i 

i=l 



Here, gi is the relative weight of the ith ob- 
servation, which is inversely proportional to a\ - 
the standard deviation -Fj(A) at wavelengths near 
the line whose profile variability is being studied. 
This definition for makes it possible to take into 
account the different qualities of the data in the 
analyzed set of spectra. TVS(X), or the Time 
Variation Spectrum [31], obeys a x 2 /{N — 1) dis- 
tribution with N — 1 degrees of freedom. 

In the case of profile variability with a suffi- 
ciently high amplitude, the TVS(X) in the vicin- 
ity of the corresponding line substantially exceeds 
its values in the adjacent continuum. To ensure 
that the increase in the amplitude is due to real 
variability, we specified a small significance level 
a <C 1 for the hypothesis that the excess is due 
to a random variation of the noise component of 
the profile. The quantites for the x 2 distribution 
are calculated in the ordinary way for this a <C 1 
level (see, for example, [32]). Let Xa be specified 



so that the probability P(x /{N - 1) > x a ) = a 
If the calculated TVS(X) value exceeds Xa, the 
hypothesis that the line profile is variable can be 
adopted. 

As was noted in [T2], for lines with small- 
amplitude profile variations or with a small num- 
ber of spectra, the function TVS(X) cannot be 
used to prove whether a given line profile is vari- 
able. In the absence of any visible variations, we 



determined whether a line profile was indeed vari- 
able using the following procedure. 

Before the standard-deviation spectrum was 
obtained, the differential spectra were pre- 
smoothed using a filter that is broad compared to 
the pixel size. If the filter width is AA on the wave- 
length scale, then the amplitude of the smoothed 
random (noise) component of the differential pro- 
files Nj(t) decreases after smoothing by a factor of 
~ \J AA/<5A, where 5X is the width of an individ- 
ual pixel on the wavelength scale in the vicinity of 
the line. If the width of the variable component is 
not smaller than that of the filter, then smoothing 
with this filter will not substantially change the 
amplitude of the variable component, while, after 
the smoothing, peaks in the standard deviation 
spectrum that correspond to the variable compo- 
nent can be detected with increased reliability. 

It was shown empirically that the best results 
are obtained for smoothing with a Gaussian fil- 
ter with a width of S = 0.5 - 1.0 A (15-30 pix.). 
To more clearly demonstrate the method used 
to search for weak profile variations, we used a 
smoothed time- variation spectrum, smTVS, which 
is a collection of spectra of the normalized stan- 
dard deviation cr(A, S) = y/TVS(X, S) with the 
differential spectra smoothed with a Gaussian fil- 
ter with a variable width S. 

Figure H] presents density diagrams for the 
time- variation spectra for the Hell A 4686 and Hg 
line profiles. Darker areas correspond to higher 
amplitudes of the smTVS spectrum. Only values 
corresponding to significance levels a < 10~ 4 for 
the hypothesis that the profiles are not variable 
are presented. The density diagrams show that 
the variability of the H^ and Hell A 4686 lines is 
clearly present. 

The smTVS spectrum for the Hell A 4686 line 
profile with the filter width S > 0.3 Aindicates 
variability of the CIII A 4673.95, Oil A 4673.75, 
4676.234, Nil A 4674.909, 4678.14, and 
Silll A 4673.273 lines, which cannot be de- 
tected using ordinary methods. With larger filter 
widths, variability of a large group of Oil, NIII, 
and Aril lines at wavelengths from. 4871 to 
4890 Abecomes visible. Note that, in spite of 
the fact that these lines are extremely weak and 
their residual intensities do not exceed the noise 
level (by one pixel along the spectrum) at the 
adjacent continuum, their variability is clearly 
detected. Unfortunately, our technique cannot be 
used to localize weak lines with variable profiles 
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Figure 3: Dynamical spectra for the profile variations for the Hell A4686 (top left), Hel A 4713 (top right), H/3 (bottom 
left), and CIIIA5696 (bottom right) lines. The interval between consecutive spectra is four minutes. Dark regions in 
the diagrams correspond to intervals brighter then the average profile (peaks), and light regions to less bright intervals 
(valleys) . 



accurately when larger filter widths are used. 

Note that the efficiency of this method of using 
the smTVS spectra to search for weak line-profile 
variability is sensitive to the number of spectra ob- 
tained, and increases substantially as this number 
increases. 

4.4 Wavelet Analysis for the Line- 
Profile Variations 

Analysis of the differential line profiles in the spec- 
tra of a star (Fig. [3]) indicates the presence of 
a number of discrete features. Many small-scale 
features are connected with the noise component 
of the profiles, whereas larger-scale features may 
be due to the regular component of the profile 



variations [12j. A very suitable mathematical tech- 
nique for studying the development of different- 
scale profile features is wavelet analysis (see, for 
example, [331 E]). It is advisable to use as the an- 
alyzing (or mother) wavelet the so- called MHAT 
wavelet, ip(x)=(L — x 2 )exp(— x 2 /2) , which displays 
a narrow energy spectrum and whose zero and first 
momenta are equal to zero. The MHAT wavelet 
is essentially the second derivative of a Gaussian 
function (taken with a minus sign) that can be 
used to describe features in the differential line 
profiles. 

Using this wavelet, the integral wavelet trans- 
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Figure 4: Time variation spectrum for the Hell A 4686 
(top) and H,3 (bottom) line profiles. The vertical axis plots 
the logarithm of the filter width in A. 

formation can be written[33, 34J: 

oo 

W(s,u) = i J f( x )^(^)dx (2) 

— oo 

oo 

= / f(x)ip su (x)dx 

— oo 

where f(x) is the studied function (in our case, 
a differential line profile). The signal energy en- 
ergy, Ew {s, u)=W 2 (s, u), characterizes the energy 
distribution of the studied signal in (s,u) = = 
(scale, coordinate) space. Further, we will ana- 
lyze the differential line profiles in the form r(V) = 
Fi(V) — Fi(y), where the Doppler shift V from the 
central frequency of the line is used as the x coordi- 
nate. In this case, the scale variable s is expressed 
in km/s. 

To study the evolution of features in the 
differential profiles, we calculated wavelet spec- 
tra, W{s,V) = W(s,V,t), for the Hell A 4686, 
HeIA4713, EL3 and the CIIIA5696 emission line 
for all times t corresponding to our spectra. We 
will call the collection of functions W(s,V,t) the 
dynamical wavelet spectrum for the differential 
line profiles. 

Figure [5] presents density diagrams for the dy- 
namical wavelet spectra of the studied lines for 
the scale parameter s = 50 km/s. The diagrams 
present the ratios d = (W — W m i n )/AW rather 
than the function W(s,V,t) itself. Here, AW = 
W max - W min , where where W max and W min are 



the maximum and minimum of W(s,V,t) for all 
possible t in the interval of V within the line pro- 
file. Darker regions correspond to larger values of 
W(s,V,t). We have marked only the values for 
the parameter d > d cu t - the cutting parameter of 
the wavelet spectrum. In Fig. [5] dcut =0.6. 

For a comparison, the average line profiles are 
placed beneath the wavelet spectra diagrams. In 
the average profiles, the dotted lines delineate the 
region — Vsini -r- Vsini for the primary, Aa 1 (see 
Table) . 

We can see in Fig. [S]that all the wavelet spec- 
tra display broad structures moving along the 
line profiles, with features moving from the vi- 
olet to red edge of the line being most promi- 
nent. The same type of features are visible in 
the wavelet spectra of the absorption Hell A 4686, 
HeIA4713 and H^ differential line profiles, mov- 
ing from -110 -r -100 km/s to « 100 110 km/s 
over the total observation time, 2 h 50 m . The es- 
timated crossing time T cross for the entire band 
of — Vsini -7- Vsini is ~ A h . Note the similarity 
between the dynamical wavelet spectrum and the 
time variation spectrum for the line profiles (see 
Section H^J. 

The HeIA4713 and Hell A 4686 lines are ba- 
sically formed in the photosphere, and their dy- 
namical wavelet spectra behave in a very similar 
way. Apart from the primary, which shifts to- 
wards the red part of the line profile, another, 
weaker, component is seen. This secondary com- 
ponent appears ~ 40 min after the start of the 
observations at ~ — 130 4- 150 km/s and moves to- 
wards more negative velocities (~ — 300 300km/s) 
by the end of the observations. This feature may 
correspond to component Ab of the 5 Ori A sys- 
tem, since profile variations associated with possi- 
ble pulsations of this star should be traced up to 
V = V m in « —300 km/s (see the Table). It is also 
possible that this feature is formed in the power- 
ful stellar wind from the primary of the 5 Ori A 
system. 

The same secondary component as that seen 
for the Hel A 4713 and Hell A 4686 lines (although 
less distinctly pronounced) is also present in the 
dynamical wavelet spectrum of the H^ line. Ap- 
preciable features of the dynamical wavelet spec- 
trum of Hg are visible even beyond the ±V sin i ~ 
300 km/s band for the most rapidly rotating star, 
Ab, probably indicating that the intense stellar 
wind from the primary contributes substantially 
to the H^ line profile. 
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Figure 5: Dynamical wavelet spectrum of the Hell A 4686 (top left), HeIA4713 (top right), YL beta (bottom left), and 
CIIIA5696 (bottom right) line-profile variations for the scale parameter 5o=50 km/s. The corresponding line profiles 
averaged over the entire set of observations are given beneath the wavelet spectra. The dotted curve indicates the region 
±V sin i for the primary component Aa 1 of the 5 Ori A system. 



The dynamical wavelet spectrum of the differ- 
ential we can see two main features, moving al- 
most parallel profiles of the CIII A 5696 emission 
line differs appreciably from those for the absorp- 
tion lines. In Fig. [5] we cam see two main features, 



moving almost parallel in the —200 -r 200 km/s 
band from the red to the violet part of the line 
profile. 

Regular structures in the wavelet spectra are 
most clearly visible in the profile of the Hel.4713 
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Figure 6: Dynamical wavelet spectrum of the HeIA4713 line profile variations for so =25 km/s and various cutting 
parameters d c ut = (a) 0, (b) 0.5, (c) 0.6, (d) 0.65, (e) 0.7, (f) 0.75. 



line. To study these structures in detail, Fig. [6] 
presents the calculated dynamical wavelet spec- 
trum for this line for the wavelet-spectrum cutting 
parameters d cu t = 0, 0.5, 0.6, 0.65, 0.7, 0,75 and 
for So = 25 km/s. As d cu t increases, more pro- 
nounced structures of the wavelet spectrum be- 
come more dominant. All the wavelet spectra dis- 



play movement of a feature from -120 to 90 km/s; 
for d cu t > 0.75, only this main feature is present 
in the wavelet spectrum. 

To clarify the nature of the features revealed in 
the wavelet spectra for the differential profiles, we 
must determine whether the line profile variations 
related to these features are regular. We consider 
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this question in the following section. 



5 DISCUSSION 



4.5 Search for Regular Variability 

To search for periodic profile variations, we carried 
out a Fourier analysis of the line-profile variability 
in the spectrum of 5 Ori A. The Fourier spectra 
were cleaned of false peaks using the CLEAN |35| 
procedure modified by Vityazev [36J. 

For each of the studied lines, we constructed 
the time series AF(t, A) of differential fluxes for a 
given profile at time t for a given A. For conve- 
nience, the A values were converted into Doppler 
shifts V from the profile center in the center-of- 
mass frame of the 5 Ori A triple system (see Sec- 
tion S3]). 

This CLEAN procedure was applied to all the 
observed V values along a line profile. To smooth 
the noise component of the profile variations, the 
AF values were averaged within a spectral win- 
dow with width AV = 5-6 km/s. Our calculations 
showed that our specific selection of AV^ do not ap- 
preciably influence the resulting Fourier spectra. 

Figure \7\ presents contour plots of the square 
of the amplitude of the Fourier transform 
(Fourier spectra) for the Hell A 4686, HeIA4713, 
CIII A 5696 and CIIIA5696 lines. When con- 
structing the density diagram for the power spec- 
trum of the Fourier transformation, all values of 
the periodogram corresponding to low significance 
levels q = 10~ 7 — 10 -3 for the hypothesis that there 
was a strong white noise peak in the periodogram 
were rejected. Thus, only values corresponding to 
significances for the hypothesis that there is a har- 
monic component in the periodogram for a value 
a = 1 — q > 0.999 are presented in the plots. 

Figure 7 has a broad peak at the frequency 
v = 5.9±l<i~ 1 . The large width of the peak is due 
to the relatively low resolution of the Fourier spec- 
trum, related to the short duration of the obser- 
vations. The value P = 1/v pa 4 ft .l ±Q h .7 exceeds 
the duration of the observations, and so cannot 
be directly identified with the period of the regu- 
lar profile variations. We will call this value the 
quasi-period for the regular profile variations. Of 
course, the presence of regular line-profile varia- 
tions with the proposed period must be confirmed 
by more extended observations. At present, we 
can only state that variations of the differential 
profiles in the spectrum of 5 Ori A are consistent 
with the period ~ 4.1±0.7 ft , with the same period 
being indicated for all the studied lines. 



Our analysis shows that regular line profile vari- 
ations with the characteristic time P = 4 h — 6 h 
may occur in the spectrum of the 5 Ori A system. 
Such periods are typical for profile variations due 
to non-radial pulsations (NRPs) in O stars [I]. 

Apart from NRPs another possible origin for 
the rapid line profile variations is rotational mod- 
ulation of the profiles. In the model of Kaper et 
al. pQ, it is assumed that, along with a spherically 
symmetrical outflow in the form of a stellar wind, 
some denser, compact corotational streams rotat- 
ing with the angular velocity of the star could be 
present in the atmosphere. It is proposed that ad- 
ditional absorption of the stellar radiation when 
these streams cross the line of sight results in reg- 
ular line profile variations. The period of these 
variations will be P n = P ro t/n, where n is the 
number of corotational streams. 

The orbital period P ro t of the primary Aa 1 
can be estimated from the data in the Table: 
-frot/sini ~ 3.6 d. The real inclination i of the 
rotation axis of the primary is probably close to 
the orbit inclination, 67 — 77° [H], which corre- 
sponds to P ro t I ' sini 3.6 d. Explaining the quasi- 
period of the variations P ~ 4^.1 as a result of 
rotational modulation of the profiles requires the 
value n = 21, which does not seem plausible. In 
addition, in this case, we should see no fewer than 
10 regular components in the entire line profile, 
which are not observed (Fig. [3]) . 

Based on this reasoning, we conclude that the 
rotation of Aa 1 cannot give rise to the detected 
rapid profile variations with the quasi-period P ~ 
4^.1. Similar estimates based on data from the 
Table indicate that the rotation of Aa 2 and Ab 
likewise cannot explain the profile variations in the 
spectrum of 5 Ori A. 

Thus, the most likely origin of the regular pro- 
file variations is photospheric NRPs of the stars in 
the 5 Ori A system. The relations I ~ 0.1 + 1.09 | 
Ac^o | /vr and m pa 1.33 + 0.54 [ Acpi \ /2pi can be 
used to determine the (1 m) mode of the NRPs, 
where A0o is the phase difference of the Fourier 
components of the line-profile variations of the 
main component z^o of the NRPs and Afa is the 
same value for its first harmonic v\ = 2v$ |30j. 
Unfortunately, due to the short duration of the 
data set A</>oand A<p\ cannot be determined with 
sufficient accuracy. The value Acpo for the rig line 
profile can be determined only in the velocity in- 
terval — 18 4- 150 km/s, and is equal to pa n. As- 
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Figure 7: Fourier power spectra for the line profile variations. Fourier spectral amplitudes are presented for low 
significance levels for the hypothesis that there is a strong peak in the white noise periodogram q ,for the Hell A 4686 
(top left; q = 10" 7 ), HeIA4713 (top right; q = 1CT 6 ), H a (bottom left; q = 10" 7 ), and CHI A 5696 (bottom right; 
q = 10 -6 ) lines. Darker regions correspond to larger amplitudes of the Fourier spectra. 



suming that this value can be extrapolated to the 
entire velocity interval ±V sin i, we find I ~ 2. 

Note that the profile variations of the studied 
lines in the spectrum of 5 Ori A are qualitatively 
similar to those in the spectra of numerous Be 
stars [38] that can be described as the result of 
NRPs in the sector mode (I, \ m |) = (2, 2). 

Let us now analyze the line profile variations 
in the spectrum of 5 Ori A supposing that they are 
related to NRPs in the (2, ±2) mode. The velocity 
of a volume element of the star undergoing NRPs 
in spherical coordinates is oc Y m i(9,(p)exp(iujt), 
where ^22(^5 <j>) (1— cos 26)*exp(im(fi) is a spher- 



ical harmonic and u> is the angular frequency of 
the pulsations [30]. The time dependence of the 
velocity is specified by the factor exp(im^ + ut), 
where 2ir/tv is the pulsation period. The phase ve- 
locity for the propagation of perturbations across 
the stellar surface is lonrp = 2vr * uj/m. This re- 
lation implies that the perturbation propagation 
period is twice the pulsation period. 

The propagation of perturbations of the ve- 
locity of matter in the photosphere corresponds 
to movement of features in the differential pro- 
files (peaks and valleys) with the same periods. 
This means that a profile feature associated with 
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NRPs will cross the band ± V sin % over for the time 
Tnrp, where T^rp is the NRP period. Therefore, 
we can write the approximate equality T^rp ps 
Across j where T cross is the time for the feature to 
cross the band ±Vs'mi in the differential pro- 
file. We found the time T cross from our analy- 
sis of the wavelet spectra of differential profiles of 
the HeIIA4686 and HeIA4713 photospheric lines 
(Fig. [5]): T cross w 3 — 4 h , very close to the quasi- 
period of the profile variations P ~ A h determined 
above. 

We conclude that the observed profile varia- 
tions are consistent with the hypothesis that they 
are associated with photospheric NRPs of the pri- 
mary, Aa , of the 5 Ori A system in the sector 
mode (I, | m |) = (2,2). Since the profile fea- 
tures in the wavelet spectra move from the violet 
to the red part of the profile, m < 0, yielding 
(l,m) = (2, -2) for the NRP mode [30]. 

Apart from these features of the dynamical 
wavelet spectra, which we suppose to be related 
to NRPs, these spectra also display weaker fea- 
tures in the violet and red line wings beyond the 
±Vsini band. These cannot be due to variations 
in the photosphere of Aa 1 . Two explanations for 
these specific features of the line profile variations 
are possible. First, they may be related to pro- 
cesses in the stellar wind in the expanding at- 
mosphere of Aa 1 . Since (i) we have estimated 
the contributions from Aa 2 and Ab to the main 
line profiles to be relatively small, (ii) CIIII A 5696 
emission appears only in the spectra of the bright- 
est hot supergiants (in our case, Aa 1 ), and (iii) the 
variability of the CIII emission is likely due to vari- 
ability of photospheric absorption lines, and also 
exceeds the limits of the ±V sin i band, we suggest 
that the variability of the main lines is primarily 
due to NRPs in the atmosphere of the Aal super- 
giant, and also to a variable contribution from the 
emission of the wind from the primary, modulated 
by its NRPs. Second, the variability beyond the 
± V sin i band could be related to the components 
Aa 2 and Ab. 

Two clearly pronounced parallel structures are 
visible in the dynamical wavelet spectrum of the 
CIII A 5696 differential line profiles (FigHJ), asso- 
ciated with the movement of broad features in 
the wavelet spectrum towards negative velocities. 
The right component moves with a velocity of 
« 180km/s to sa 40km/s and is located fully 
(taking into account its own line-profile widths) 
within the ±Vsinz band. The band crossing 



time for this component is ~ 6 , which substan- 
tially exceeds the crossing time for the compo- 
nents of the Hell A 4686 and HeIA4713 lines re- 
lated to NRPs. These structures are undoubtedly 
associated with some imhomogeneities in the at- 
mosphere that contribute to the emission at the 
CIII A 5696 line frequencies, rather than directly 
with NRPs of the primary. The possibility of gen- 
erating quasi-regular structures in a stellar wind 
originating from photospheric NRPs is indicated 
in [37]. 

In support of this idea, we note that the pat- 
terns of the dynamical wavelet spectra related 
to gaseous imhomogeneities detected in the CIII 
emission, and also the wavelet spectrum features 
seen in other lines beyond the ± Vsini band, are 
mutually complementary. Bringing into coinci- 
dence the wavelet spectra for the HeIA4713 pho- 
tospheric line and the CIII A 5696 envelope line 
(Fig. [5]) indicates an avoidance of the wavelet 
spectra; i.e. the absence of any features in the 
CIIIA5696 emission wavelet spectrum in the veloc- 
ity range containing the features in the HeIA4713 
photospheric absorption wavelet spectrum. 

To clarify the extent to which the obtained 
quasi-period for the line-profile variations in the 
spectrum of 5 Ori A corresponds to possible pul- 
sation periods of OB stars, we have plotted the 
positions of Aa 1 and Aa 1 on an diagram of ef- 
fective temperature vs. pulsation period in the 
quadrupole mode I = 2 for hot main-sequence 
stars ( [39J, Figs. 6 and 7]. The positions of 
the components in Fig. [8] (top) are below the 
short-period boundary of the pulsation-instability 
strip for these stars, and do not correspond to the 
masses of the stars determined in[T3] (see Table). 
At the same time, the position of Aal in the di- 
agram calculated for a high-mass main-sequence 
star with M =30 M. (Fig. [HI bottom) accurately 
corresponds to the pulsation-instability strip for 
this star. 

We suggest that, even in the case of substan- 
tial mass- loss from a star, its inner structure does 
not vary appreciably, and the periods of its NRPs 
correspond to pulsations for a main-sequence star 
whose mass is equal to the initial mass of the star 
up to the mass-loss stage. 

Confirmation of the reality of the obtained pe- 
riod and clarification of the nature of the line- 
profile variations requires spectra of 6 Ori A ob- 
tained over two to three nights, making it possible 
to encompass four to six NRP cycles. 
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Figure 8: Position of components Aa 1 and Aa 2 (grey circles) in a diagram of logT e g vs. pulsation period P [35]. The 
top plot shows the periods for the pulsation modes I = 2 for (3 Cep type stars and slowly pulsating B stars; the small 
dots indicate the calculated periods, and the bold dots the pulsation periods of main-sequence stars with masses of 5Mq 
and 12Mq . The bottom plot shows the periods for pulsation modes / = 2 for a main-sequence star with a mass of 30Mq . 



6 CONCLUSION 

The following conclusions may be drawn from our 
observations. 

1. All the studied lines display variable pro- 
files, with the variability amplitude being 0.5%- 
1%. This conclusion has validity level of 0.999. 

2. Large-scale components in the velocity in- 
terval 25-50 km/s are detected in the dynamical 
wavelet spectrum of the Hell A4686, Hel A4713,H., 
and CIII A5696 line profile variations, which move 
in the — V sin i-r- V sin i band for the primary of the 
system, Aa 1 , with a band crossing time of 4-5 h . 
Some variable features that go beyond the band 
are most likely due to emission imhomogeneities in 
the wind from SOvi Aa 1 . The less bright compo- 
nents of the system Aa 2 or Ab may also contribute 
to the variability beyond the — Vsini 4- Vsini 
band. 3. We have found short-period variations 
of the studied line profiles with the characteris- 
tic times ph 4 h .l.. We present evidence suggest- 
ing that these variations are associated with non- 
radial pulsations of the 09.511 primary of the sys- 
tem, Aa 1 , in the sector mode (l,m) = (2, —2). 
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